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Abstract 

Lipid bilayer membranes below their main transition have two tilt order 
parameters, corresponding to the two monolayers. These two tilts may be 
strongly coupled to membrane shape but only weakly coupled to each other. 
We discuss some implications of this observation for rippled and saddle phases, 
bilayer tubules, and bicontinuous phases. Tilt difference introduces a length 
scale into the elastic theory of tilted fluid membranes. It can drive an instabil- 
ity of the flat phase; it also provides a simple mechanism for the spontaneous 
breaking of inversion symmetry seen in some recent experiments. 

The curvature model of fluid bilayer membranes has proved quite successful in explaining 
the shapes of membranes above their main transition In this model locality, coordinate 
invariance, and bilayer symmetry restrict the form of the energy functional for shapes to 
just two terms, involving the mean and Gauss curvature. In fixed topology the total Gauss 
curvature is constant, and so the minimum-energy conformation is a surface of vanishing 
mean curvature, for example a fiat plane. 



Below the main transition, additional degrees of freedom enter the elasticity of mem- 
branes as their hydrocarbon chains begin to order. In analogy to smectic liquid crystals, 
one expects a soft tilt degree of freedom to appear, reflecting the spontaneous breaking of 
rotational invariance in the plane. Helfrich and Prost began the systematic study of the 
mutual influence of tilt order and membrane shape 0. A number of nonflat ground state 
phenomena found in membranes below their main transition have since been attributed to 
tilt, including rippled phases and tubule phases (see for example [IHl])- Tilt order also 
proves to be crucial for the intrinsic chimlity of individual amphiphiles (if any) to influence 
the conformations of membranes 0. 

Despite much progress, however, a number of mysteries remain in the study of one- 
component, symmetric bilayer membrane conformations. For example, experiments with 
racemic mixtures of lipids have found asymmetric ripple ground states, which are chiral . 
Similarly, achiral amphiphiles can form tubules , which again appears to require chirality 
Even chiral lipids easily form helical ribbons of either handedness; sometimes a single 
ribbon appears to switch handedness in the middle of its growth |]TU[. Finally, cubic phases 



of bilayer membrane are predicted to be scale-invariant in the pure curvature model 
One might therefore expect them to collapse to a microscopic cell size, but in fact they can 
be stabilized at well-defined, mesoscopic scales . Theoretically, selection of a length-scale 



of several nanometers has been attributed to higher order curvature terms both for cubic 
phases |jTl| and for a presumed hats and saddle super-structure of fluid bilayers |jl3 . 



In this letter we explore a new model for the conformation of membranes: we augment the 
curvature model with two tilt director fields corresponding to the two monolayers. Thus our 
work fits into the general program of taking more seriously the bilayer aspect of membranes, 
both in their equilibria and dynamics (e.g. p^-p!6|). When the directors are aligned we 
reproduce existing models; when they are not we get new physics. Imagining that the tilt in 
a monolayer induces spontaneous curvature orthogonal to its direction, in the corresponding 
bilayer the flat state with anti-parallel tilt then is frustrated and thus more costly than 
a saddle conformation in which the tilts are oriented orthogonal to each other. Without 



interaction between the two layers the flat state is always unstable locally against this kind 
of saddle conformation. Real membranes will always have at least some tendency to align 
the tilts. The instability will then occur only if the anisotropic spontaneous curvature is 
larger than a threshold value determined below. 

There seems to be no direct experimental evidence for or against non-aligned tilt in 
bilayer phases. We will motivate the model, show how it gives a new mechanism for the 
transition from fiat to rippled membranes, and comment on how it addresses the questions 
above. 

Our work was first motivated by a desire to understand the origin of a nonanalytic 



curvature energy proposed by Fischer [0. Our model is mathematically similar to one 



independently proposed by Fournier [|T8|, but the physical motivation is quite different: while 
he considered an anisotropic impurity adsorbed onto a membrane, our tilt is an intrinsic 
property of a pure bilayer, and hence quite generic. Other differences will be noted below. 

Model: We will restrict attention to systems of nonchiral amphiphiles. In the covariant 
notation developed in 0,|T^,^ , this means we consider only elastic energy terms constructed 
without the in-plane antisjTumetric tensor eap- To focus attention on the new elements we 
will also impose an additional "nematic" symmetry (see below). This assumption is strictly 
for mathematical simplicity; we leave the full model to future work. 

Above the main transition we imagine the membrane to be two identical 2-dimensional 
fluid sheets of elastic monomers, independent except for the constraint that they lie a fixed 
distance above and below a common surface. We will label the layers arbitrarily as "+" and 
"— " , but since the layers are identical we will insist that nothing changes if we reverse the 
labeling. To define the curvature tensor Kajj we will choose the normal vector N pointing 
from the "— " to the "+" side. Each sheet has its own bending stiffness, stretching modu- 
lus, and spontaneous curvature. When we combine the sheets, the bending stiffnesses and 
stretching moduli add, while the spontaneous curvatures cancel (for details see e.g. [0). 
We will neglect the stretchiness of the membranes, leaving only the curvature stiffness. 
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Below the main transition each monolayer develops a local average tilt. We will not be 
interested in the main transition itself, so we will take the director of the amphiphiles to be 
at a fixed angle to the layer normal, and the degree of ordering to be constant. In other 
words, we will describe the tilt by a pair of unit vector fields m± tangent to the surface: 
these are the normalized projections of the average molecular directors. The elastic energy 
is then a local functional of m+, m_, and the membrane shape, described by its curvature 
tensor Kajj- As mentioned above, we will rather artificially assume invariance when either 
m+ or m_ changes sign. 

Traditionally one takes m_|_ = — m_ (e.g. see |jl9|), or in other words assumes that the 
average directors in each layer are coUinear. The reasoning is that while an overall rotation 
of both m-i- is related to a broken symmetry, still the relative angle is not, and so is expected 
to lock to a preferred value. We propose to explore what happens when this assumption 
is relaxed. Our motivation is the observation that the degree of interdigitation of the lipid 
chains between the layers is in fact quite small, as deduced for example from measurements 



of the interlayer friction coefficient WM. We will allow for some weak aligning potential 
g{{m^ ■ m_)^), but as we will show other energetic contributions can readily overcome it. 
For simplicity we will continue to assume that g is minimized at = — m_, but later on 



we will reconsider this assumption too. In ref. ||T8[, g was taken to be zero, an extreme limit 
of our model. 

While the direct interaction of the tilts may be small, there will certainly be an important 
indirect interaction mediated by the coupling of tilt to the common shape. To lowest order in 
the curvature tensor the effect of tilt on a monolayer is to create an anisotropic spontaneous 
curvature, by adding the terms 

/± = /? K^^mlmi = /3 m± ■ K ■ m± (1) 



to the elastic free energy. (3 is a new parameter depending on the degree of ordering . 
The combination — is then invariant under renaming the two leaves of the bilayer. It 
vanishes if m_|_ = — m_, but more generally we need to keep these terms. At higher order 



in curvature we have the more famihar terms K^{Kcipm'^m±Y + KiKa/sK/s^m'^m']- (see [^ ). 
Since (||) will already drive our instability we will neglect the higher-curvature ^3, K4 terms. 
All told, then, our simplified model is defined by the elastic energy density functional 

/ = f + fK^pK^P + (?((m+ ■ m_)2) 

+ [±/3i^a/3m"m^ + ^Vam^V"m/3 + f (V«m")2] . (2) 

m=m-|- 

The constants ni, K2 are related to the usual mean and Gaussian rigidities, while ki, /c2 are 
related to the usual rigidities of the XY model. 

A remarkable feature of (0) is that the parameter [3 has dimensions of inverse length. 
This is a key qualitative difference from the case of parallel tilt directors: achiral symmetric 
bilayers with one tilt admit only dimensionless couplings [^,0 . One might expect that the 
value of such a parameter would be very large, on the order of the microscopic wavenumber, 
and would create correspondingly small structures. From its geometric origin, however, it 
is easy to see how (3 could be smaller than naive dimensional analysis would suggest. For 
example, the average tilt could be rather small, leading to a small net spontaneous curvature. 
In this way our term could set a mesoscopic scale. 

Stability of Bat surfaces: To understand the physics of the model (0), we first note that if 
the director in the "+" layer is pointing in the x direction, and the one in the "— " layer in the 
y direction, the (3 term becomes f3[Kxx — Kyy], which favors saddles (or other non-spherical 
shapes) . 

For a quantitative stability analysis of the fiat state in which both nematic fields are 
parallel, we need an explicit form of the interaction g. A simple form that favors parallel 
alignment is g = 7(1 — (m+ ■ m_)^)/2. We parametrize small deviations from this parallel 
state for the director fields as m± = (j!"!'/^) ~ ±0^^)- ^^^^ representation, the coupling 
term between the two layers becomes g = 270^. 

Adding up all energies to quadratic order in a Fourier representation for (j){x, y) and the 
height h{x,y) leads to 
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(3) 



with K = Ki + K2. Minimizing with respect to /iq yields /iq 
result into (0) leads to the effective energy for the (p field as 




0q. Inserting this 



/ 




2 



(4) 



The stability criterion for the fiat phase is now obvious. For weak enough interlayer 
coupling, 



the fiat phase becomes unstable to a long- wavelength modulation. This is our main result. 
The preferred directions for the q vectors are Qx = ^Qy = ±v^q^/2. Adding just two modes 
with qi = — q2 yields a ripple shape whereas adding four modes with = ±qy leads to a 
egg-carton-like square modulated phase. 

Lubensky and MacKintosh also obtained symmetric ripples in a nonchiral model. Then 
balanced an effectively negative (Vm)^ term against a stabilizing (V^m)^ term ^j. Far from 
the main transition such a balance is likely to select a microscopic length scale. In contrast, 
we have seen how our model can select long lengths. Other models assumed hexatic order 
IPJ^, while we have not. 

Beyond instability: In order to distinguish the two alternatives, ripple or saddle, we have 
to go beyond the instability. We will do so in the following using exact analysis, a simple 
variational shape and numerical minimization. To keep our formulas tractable, in this section 
we will neglect the tilt stiffness terms, i.e. we set ki = k2 = 0. Thus the tilt fields track 
the curvature exactly, and in particular can change discontinuously by 7r/2 when the mean 
curvature changes sign. A more realistic model would broaden these discontinuities into 
domain walls, with an energy cost per unit length. 

First, we discuss the ripple phase. For a one-dimensionally modulated conformation, the 
free energy density takes the form 



7 < PV^, 



(5) 
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/ = (k/2) + /3/i^^(cos2 0+ - cos2 (/,-) 

+ ^(1 — (cos0^ COS (p^ + sin 0+ sin 0^)^). (6) 

One could try as a trial variational function a sinusoidal ripple with wavenumber q and 
height /ig and nematic fields always arranged perpendicular to each other which leads to 
an energy density of f/n — —^{P'^^ ~ l)- However one can easily see that the sinusoidal 
shape is far from optimal. If the tilts are orthogonal, thus rendering the /5 term —\(3hxx\, 
ripples that consist of circular arcs of radius k/P are lower in energy than sinusoidal ripples. 
A straightforward calculation shows that allowing the fields to assume a non-orthogonal 
configuration docs not lead to a lower energy for 7 < P'^/k, which is the region of interest. 
Thus, in summary, as long as only one-dimensionally modulated phases are considered, a 
ripple phase consisting of circular arcs with the nematic fields orthogonal in the two layers 
is the most favorable configuration. Due to the absence of gradient terms for the director 
fields, these fields jump whenever the curvature changes sign. Including the gradient terms 
would thicken the jump region to a domain wall; the energy of this wall would then favor 
long wavelength (fewer walls per unit length), breaking the degeneracy we have found so 
far. 

When modulation of the membrane in two dimensions is allowed, the free energy density 
takes the form 

/ = {k/2){K^ + hyyf + (7/2)(l - cos2(0+ - 0-)) 
+ P[hxx{cos^ 0+ — cos^ + hyy{sin'^ 0+ — sin^ 0~) 

+ /ia;2/(sin20+ - sin20-)]. (7) 

Again a variational trial function sinusoidal in two directions is not the best choice. It leads 
to an energy density oi f / k = — — 7)- Using the fact that the lowest energy ripple 

conformation had constant curvature, we construct a saddle phase in which the curvature 
is piecewise constant. This leads to a saddle conformation in which the corrugations are 
parabolic arcs. In regions where the principal curvatures have the same sign, we allow the 
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nematic fields to be parallel (thus gaining from the 7 term), and where they have opposite 
signs the fields remain orthogonal. The free energy density of this state is f / k = — — 7). 
This energy is never lower than either the fiat phase or the circular arc ripples even though 
it improves the value of the sinusoidal saddle for a large range of 7 values. 

Having exhausted the range of analytic saddle conformations, we turned to a numerical 
minimization of the full free energy (^. While we could improve on the parabolic saddles, 
still we found no saddle phase with lower energy than either the circular ripples or the flat 
phase. 

Discussion: We have not obtained a unique wave-length for the instability. What we found 
was rather that the tilt-difference coefficient f3 sets a preferred curvature for cylindrical 
segments. For small amplitudes this means that the combination hgq'^, but not q itself, is 
fixed. This happened because we neglected the gradient terms in the tilt fields: the energy of 
a configuration is then invariant under a rescaling of h{x, y) —>■ Xh{Xx, Xy) if at the same time 
7 — i> 7/A^ and (3 P/X. Gradient terms in m, of course, would favor large wave-length, up 
to a maximum of ~ k//?, where the membrane rolls up into cylinders. One way to stabilize 
finite-wavelength ripples would be to stretch the membrane with a lateral tension. 

Alternatively we can take our results as indicating an instability to a phase of cylinders 
or other curved objects, especially saddles. Particularly intriguing is the possibility that the 
tilt-difference term could set a scale for bicontinuous phases without resorting to microscopic 
(higher-curvature) terms as in [jll] . 



We now consider briefly the possibility of breaking spatial inversion symmetry. Other 
authors have obtained rippled or striped phases in tilted bilayers by explicitly breaking parity 
invariance (e.g. P,^). In particular, explicit parity- violating terms in the free energy have 
so far been necessary to obtain asymmetric, or sawtooth, rippled phases. But the presence 
of chiral amphiphiles does not guarantee a parity-breaking Landau energy |2^. Moreover 
as mentioned earlier, several arguments point to the possibility of spontaneous breaking of 
parity invariance: the Landau energy has parity invariance, but its minima do not. 
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In an elegant paper, Selinger et al. have proposed possible mechanisms for spontaneous 
parity breaking in monolayers Two of these could also apply to pure bilayers: (a) if 
the membrane has both tilt and hexatic order, the corresponding directors may lock to a 
fixed, nonzero relative angle; (b) conceivably two distinct local packings of molecules could 
be preferred, each of which is the other's mirror image. We would like to point out a very 
simple, concrete option with some elements of each of these: (c) a tilted membrane may 
prefer on packing grounds to order its two tilt directors at a fixed relative angle a = ±ao, 
i.e. m_|_ ■ m_ = cosoq. This option may prove more generic than the ones above. 

To see how our proposal leads to parity breaking, note that ip = (m_|_ x m_) • N is a 

well-defined pseudoscalar order parameter, ip vanishes as a — 0, as required. 

cos a/2 ± sin a/2 

=l=sina;/2 cosa;/2 
1 
-1 



We can write m-i- 



Writing 1 for the unit matrix and e 



m, where m is a common tilt variable. 



the terms (|l]) now read 



/+ - = /?m • [(cos f 1 - sin f e)K(cos f 1 + sin f 
— (cos f 1 + sin |£)k(cos |1 — sin ■ m . 



(8) 



Rewriting this as (2/3 sin ^)m- [eK] -m, we recognize the chiral tilt-shape coupling introduced 
by Helfrich and Prost 0, in the coordinate- invariant form given in ||19|. This is precisely 
the term responsible for the formation of asymmetric ripples in the work of Lubensky and 
MacKintosh ||^. A similar analysis, dropping the requirement of nematic symmetry, could 
also yield the chiral term which gives rise to tubules and helices [^. 

In conclusion, we have shown that tilt difference may be expected quite generally to 
affect the conformations of symmetric bilayer membranes, introducing a new intermediate 
length scale and favoring cylindrical and saddle curvature over flat or spherical shapes. Tilt 
difference could also provide an attractively general mechanism for the spontaneous breaking 
of parity invariance. 
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